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ON GAUSS-LOBATTO INTEGRATION ON THE TRIANGLE 



YUAN XU 



^_^ Abstract. A recent result in [2] on the non-existence of Gauss-Lobatto cu- 

y—^ bature rules on the triangle is strengthened by establishing a lower bound for 

f^ the number of nodes of such rules. A method of constructing Lobatto type 

^sj cubature rules on the triangle is given and used to construct several examples. 



1. Introduction 



Recently in |6], motivated by hp-faiite element method, Gauss-Lobatto cubature 
rule on the triangle 
1^ A:^{{x,y):0<x,y,x + y<l} 

is studied, which requires {n + l){n + 2)/2 nodes, with {n— l)(n — 2)/2 nodes in the 

interior, n—1 nodes in each side and 1 point at each vertex, and is capable of exactly 

C^ integrating polynomials of degree 2n — 1. We call such a rule strict Gauss-Lobatto. 

r^ The main result of [6 shows that such rules do not exist. In this note we consider 

^^ cubature rules that have nodes in both interior and boundary of the triangle, and 

establish a lower bound for the number of nodes, from which the non-existence 

J> of the strict Gauss-Lobatto rule follows immediately. We also study the structure 

00 of rules that attain our lower bound and give a method for constructing cubature 

'^ rules of degree 2n — 1 with n—1 nodes on each side and 1 node at each vertex. 

~^ The development is based on the observation that cubature rules with nodes on 

. the boundary can be constructed, by restricting to the class of bubble functions 

[^ (functions that vanish on the boundary), from cubature rules with nodes on the 

interior. This leads to a bootstrapping scheme for transforming some cubature rules 
with interior points into higher order rules with a specific number of nodes on the 
boundary. Several examples are constructed to illustrate the algorithm. 



k> 2. Results 

^ For n G No, let II^j denote the space of polynomials of (total) degree n in two 

variables. It is known that 

,. tt2 fn + 2\ (n + l)(n + 2) 
dimn„ - ' ' - ^ 



2 y 2 

Let W{x,y) be a non-negative weight function on the triangle A with finite mo- 
ments. A cubature rule of precision s with respect to ly is a finite sum that satisfies 

N 

(2.1) / f{x,y)W{x,y)dxdy^y2Xkf{xk,yk), V/ e R^. 

^^ fc=i 
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We choose W to be the Jacobi weight Wa,i3,'j{x, y) — x°'y^{l ~x~ y^ for a,(3,j> 
— 1. The case a — (3 — j = corresponds to the constant weight. These weight 
functions are often considered together with the orthogonal polynomials, called 
Jacobi polynomials on the triangle, that are orthogonal with respect to them; see, 
for example, [3 p. 86], and [1, 2 in connection with cubature rules. 



For the Jacobi weight, it is known that the number of nodes N for (2.1 1 satisfies 



(-) N>_m-^ii ■'---' 



2 ' L2- 

if s^2n- 2. 



This lower bound is classical for s = 2n — 2 ([8]) and given in [T] for s = 2n — 1, 
which agrees with Moller's lower bound for centrally symmetric weight functions [7\ 
as well as [5]. A cubature rule that attains the lower bound is naturally minimal, 
meaning that it has the smallest number of nodes among all cubature rules of the 
same degree. The lower bound, however, is most likely not sharp; that is, a minimal 
cubature could require more points than what the lower bound indicates. Minimal 
cubature rules are sometimes called Gaussian cubature rules. Their construction is 
closely related to orthogonal polynomials of several variables. For discussion along 
this line, see [31 [3] and references therein. 

The cubature rules that we consider are of precision s and are of the form 

„ No Ni 

(2.3) / f{x,y)Wa,0,'fix,y)dxdy = ^Xk.of{xk,n,yk.o) + ^>'k,ifixk.i,0) 
■'^ k=i k=i 

N2 N3 

+ X! ^fc,2/(0, yk.2) + X! ^k,3f{xk,3, 1 - Xk,3) 
k=l k=l 

+ Mo/(0,0)+Aii/(1,0)+M2/(0,1), 

where {xk.o,yk,o) are distinct points in the interior of A, (a;fc_i,0), (0, j/fe,2), and 
(xfc^s, 1 — Xk,3) are distinct points on the side y = 0, x — 0, and x + y = 1 (but not 
on the corners) of A, respectively, Xk.j > and fii > 0. Such a cubature has 

N ■.^No + Ni+N2 + N3 + 3 

nodes. The main result in ^ states that such a cubature rule does not exist if 
s = 2n — 1 and 

which has a total number of nodes N = {n + 2){n + l)/2. This follows as an 
immediate corollary of the following theorem. 



Theorem 2.1. If a cubature rule of the form ( 2.3 1 exists with precision s = 2n— 1 
or s — 2n, then 



(2.4) No > 



i^ ifs = 2n-l, 



n{n— 1) 



2 ' L^J »/s = 2n-l. 



(2.5) iVo + iV.>|^ -^ ;^^^^^^ ' .^1,2,3. 
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Proof. The cubature (2.3) exactly integrates degree s polynomials of the form 



(2.6) 



xy{l-x-y)g{x,y) if 

g{x,y)Wa+i,i3+i,j+i{x,y)dxdy 



A 



No 

E 
fc=i 



K,o9ixk,o,yk.Q), V.gen^ 



-3' 



where A^ q = ^k.oXkfiVk.oi^ — Xk.o — yk.o), which is a cubature rule of precision 
s — 3 for the weight function Wa+i./s+i.-y+i so that, by (2.2), iVo has to satisfy 



the lower bound in the inequality of (2.4). On the other hand, the cubature (2.3) 



exactly integrates degree s polynomials of the form f{x, y) = x{l — x — y)g{x, y) if 

V5en2 



-3j 



A 



g{x, y)Wa+i.j3,'i+i{x, y)dxdy 



No 

y^,^k,Qg{xk,o,yk,o) 
fe=i 



Ni 

E 

k=l 



KA9ixk,i,0) 



where Afc,o = >^kfiXkA^ - ^k,o - Vk,o) and A^ ^^ = Afe,ia;fc,i(l - Xk,i), which is a 
cubature rule of precision s — 2 for the weight function Wa+i,i3,-y+i{x,y)^ so that 
A'o + iVi satisfies the lower bound in (2.2), which gives the inequality of (2.5) for 
i = 1. Similarly, we can derive lower bound for A'o + N2 and Ni^ + N^. D 

One naturally asks if there is any cubature rule that attains the lower bound in 
the theorem. For s = 2n — 1, this asks if there is a cubature rule of precision 2n — 1 
with 



(2.7) 



N^ = 



n{n — 1) 



and Ni = 



n — 1 



i = l,2,3. 



We expect that the answer is negative. A heuristic argument can be given as follows: 
Assume that A'o — "^"„~ ' . Then the proof of the theorem shows that (2.6) is a 



cubature of degree 2n — 4 with n{n — l)/2 nodes, which is known to exist only for 
small n. Assume that it does exist. We define a linear functional £1, acting on 
polynomials of one variable, by 



(2.8) Cig 



g{x)W, 



A 



a+1,/3,7+1 



'x,y)dxdy- 



No 

E 
fe=i 



Afe,oa;fe,o(l - Xkfi - yk,o)gixk,o), 



where Xk,0: 2/fe,o and Afc,o are as in (2.3 1. Applying ( |2.3[ ) on polynomials of the form 
f{x, y) = g{x)x{l — X — y) shows that 



(2.9) 



Nt 



Vg e n2„- 



where A^ j^ = Xk,iXk,i{^~Xk,i)- The functional Ci defined in equation (2.8) defines 



a bilinear form \p,q\i = Ci{pq) that could be indefinite ([(?,g]i is not necessarily 
positive). If the bilinear form were positive definite on n2„_3, then (2.9) could 
be regarded as a quadrature rule of A'l nodes and of degree 2n — 3 for £1 and, 
consequently, A'l > n — 1 by the standard result in Gaussian quadrature rule. 



which is stronger than the second equation of (2.7). Thus, in order for (2.7) to 
hold, we would need £1 to be indefinite on n2„-3, that is, C{q^) = for some 
nonzero q £ n„_i, and we would need to require that £1 has a quadrature rule 
of degree 2n — 3 with A'l = [^^y^J nodes. A simple count of variables (the nodes 
and weights of the quadrature rule) and restraints (the polynomials that need to 
be exactly integrated) shows that this is unlikely to happen, although it still might 
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as the equations are nonlinear. For a linear functional that defines an indefinite 
bilinear form, the theory of Gaussian quadrature rule breaks down since orthogonal 
polynomials may not exist and, even they do, they may not have real or simple 
zeros. In particular, we do not have a lower bound for the number of nodes of a 
quadrature rule for such a linear functional. 



The above argument indicates that if we want the cubature rule (2.31 to have 



the smallest number of interior points, then the best that we can hope for will be 
a cubature rule of degree 2n — 1 that satisfies 

n{n — 1) 



(2.10) 



Nn = 



and N,>n~l, I < i < 3. 



We shall call a cubature rule that attains the lower bound in (2.10) Gauss-Lobatto 



cubature rule. The proof of the lower bound indicates how such a cubature rule can 
be constructed. Let us also define linear functionals 



No 



^29 ■■= / g{y)Wa,i3+i,^+i{x,y)dxdy ~ > Afc,o2//c,o(l - Xk.a - yk,Q)g{xkfi) 



A 



^39 



g{x)W, 



A 



a+1,/3+1,7 



{x,y)dxdy- 



fc=i 

No 

E 
fc=i 



^k,oXkflykfi9{xk,o) 



Then we can summarize the method of construction as follows. 

Algorithm. We can follow the following procedure to construct a Gauss-Lobatto 



cubature rule of form (2.3) 



Step 1. Construct a cubature rule of degree s — 3 for Wa+i.p+i.-y+i in the form 
of (2.6) with all nodes in the interior of A, and define 

A* 
(2.11) Afe,o 



^*;,0 



(a;fc,o2/fc,o(l -a;fe,o-yfe.o))' 



Step 2. Construct Gaussian quadrature rules (2.9) and 

N2 Ns 

^29 ^ Yl ^*k,29{yk,2), C39 == Y^ K,39{xk,3), V5 e n2„_3, 

fc=l k=l 

with respect to the linear functional Ci, £2,^3, which gives the nodes of the cuba- 



ture rule (2.3) on the boundary of the triangle, and define 
(2.12) Xk^i = - — ^_, — ^, Xk,2 = - — /-, — 7, Afc^3 



A 



A;,3 



Xk,i{'^ - Xk,i)' ■' yfe,2(l-2/fe,2)' ^ a;fc^3(l-a;fe_3)' 

Step 3. Finally, the weight /io, /ii, and ^2 are determined by setting f{x,y) — 
l,x,y in (2.3) and solve the resulted linear system of equations. D 



Proposition 2.2. Assume, in the above algorithm, that the linear functional Ci, C2 
and C3 are positive definite on H2n~3 o-nd all Xk,i,yk,2,Xk^3 are inside (0, 1). Then 
the algorithm produces a cubature rule of degree 2n — 1 in the form (2.3). In 



particular, if Nq — 



then the cubature rule is Gauss-Lobatto. 



Proof. We need to show that the cubature rule constructed by the algorithm holds 
for aU / e Tlln-iy that is, ( p^] ) holds for all / e LE^ with s = 2n - 1. A moment of 
reflection shows that, with z = 1 — x — y,lij. can be decomposed into a direct sum. 



(2.13) 



Ul = xyzUl 



xzU,_2[x] © yzUs-2[y] ® xyn,_2[x] ® Ul, 
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where ns_2[x] and ns_2[2/] denote the space of polynomials of one variable in x 
variable and y-variable, respectively. If / G xyzn^__3, then (2.3) reduces to (2.61, 



which holds by our construction. If / € xzWs-2[x\, then (2.3) reduces to (2.9) 



for £i, which holds by our construction. The same holds for / e yzY\s-2[y\ and 
/ S xyIis-2[x\^ whereas for / e 11^, the cubature is verified by Step 3. Thus, by 
(2.13), the cubature rule holds for all / € n2„-i. D 



It should be pointed out that, as long as we have a cubature rule of degree 2n — 4 
with all nodes in the interior of A in the Step 1, regardless if it is a minimal one, 
the Step 2 and Step 3 could be carried out and Proposition |2.2| applies. Thus, the 



algorithm can be used to construct cubature rules of degree 2n—l in the form ( 2.3 ) 
with Ni= N2^ N^^n-1. 

Let us comment on how the steps in the algorithm can be realized. 

For Step 1, a cubature with the specification can be constructed by solving 
moment equations, that is, solving the system of equations formed by setting 



g{x,y) — x'^y^ for i+ j < s — 3 in (2.6) for Xkfi, yk,o and A^g. There have been 
a number of papers based on this method, see e.g. [TT] for the latest result and 
further references. Another approach is to use a characterization of the cubature 



rules that attain lower bound in (2.4), which is given in terms of common zeros 



of certain orthogonal polynomials and can be used to find cubature rules of lower 
order, see [71 H]. Not all cubature rules obtained via either methods work for our 
purpose, since we require that all nodes are inside the domain. 

For Step 2, one needs to check that Ci,C2,C3 are positive linear functionals on 
II2TI-3- Once they are, the standard procedure of constructing Gaussian quadrature 
rules applies. In particular, we can apply the standard algorithm to generate a 
sequence of orthogonal polynomials up to degree n—1 with respect to Ci inductively; 
the nodes of the Guassian quadrature rule of degree 2n — 3 for Ci are the zeros of 
the orthogonal polynomial of degree n — 1 with respect to Ci . 

In the following we give several examples of Gauss-Lobatto cubature rules of 
degree 5 and 7 for the unit weight function. The minimal cubature rules for the 
degree 5 and 7 have nodes 7 and 12, respectively (f?^). For Gauss-Lobatto rules. 



the number of nodes are necessarily larger, as seen in (2.7). 



Example 1. There is a cubature formula of degree 5 with 12 nodes, 3 interior, 2 
on each side and 1 at each corner of the triangle. To illustrate our procedure, we 
shall present the nodes and weights in steps. 

The three interior points and weights are given in the first table. While the nodes 
are those of a cubature rule of degree 2 for the weight function M^i,i,i(x, y) = xy{l — 
X — y), the weights are relates to those of the latter cubature rule by (2.11 ). These 
nodes are common zeros of quasi-orthogonal polynomials (see |10| for definition) of 
degree 2 which were found by solving the nonlinear system of equations in Theorem 
4.1 of HOl. 



Xk,0 


Vkfi 


Afe,0 


0.15881702219143 


0.19201873632215 


0.101342396527698 


0.56219234596964 


0.19201873632215 


0.117181247909596 


0.22100936816107 


0.55798126367785 


0.118066904793533 
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where the nodes are the common zeros of foUowing three polynomials of degree 2: 
19 - Vl05 + x(-91 + Vl05 + 112a;) + 2(-7 + %/l05)y, 
49 - 3\/l05 - 175a; + 5\/l05a; + 112a;^ - 84?/ + 4Vl05y + 224:xy, 
154 - 6\/l05 - 301a; + ll%/l05a; + 112a;2 - 609y + 7%/l05y + 560a;?; + 560?/^ 

The nodes on the edges of the triangle, but not on the corners, are the nodes of 
Gaussian quadrature rules for £i, £2, 'Ca, respectively, where the weights are related 



to those of Gaussian quadrature rules by (2.12). They are given in the following 
three tables: 



Xk,l 


Afe.i 


0.3931870086016 


0.02991955921794 


0.8595419130359 


0.01756588222187 



where the nodes are zeros of the orthogonal polynomials pi of degree 2 for £1 



Pi{x) 



448 V 



469 - 9V105 



)"+44k( 



889 + 61V105 



yk,2 


Afc,2 


0.4305843026985 


0.02290932968619 


0.7924406473476 


0.02022650113138 



where the nodes are zeros of the orthogonal polynomials p2 of degree 2 for £2, 
P2{x) = ^^ + ^ (-29 + VT05) a; + ^ (63 + %/l05) ; 



Xk,3 


Afc,3 


0.2629899118578 


0.02514330117112 


0.7030163143652 


0.03109870484395 



where the nodes are zeros of the orthogonal polynomials ps of degree 2 for £3 



^3(2;) 



3843 V 



10997 + 51V105 X 



10843 
Finally, the weights /ii, /i2, /i3 are given by 



)98 V 



3098 



665 - 9V105 



Ml 

0.0081170837035 



M2 

0.00326155091683 



M3 

0.00516753787639 



D 

The nodes of this Lobatto cubature rule are depicted in Figure 1. 
Our next two examples are Gauss-Lobatto cubature rules that are symmetric 
in the sense that the nodes are invariant under the symmetric group of A, or the 
permutation of {x,y,l — x — y). A symmetric Gauss-Lobatto rule for the constant 
weight takes the form 

„ Ma 

(2-14) / /(a;, y)dxdy = Y] Ak [fiuk,o, Vk.o) + /(«/c,o, Wk,o) + fiwk.o, "fe.o)] 

Ml 

+ Y,Bk [fiuk,i,0) + /(0, 1 - Uk,i) + /(I - Uk,i,Uk)] 

k=l 

+ C[/(0,0) + /(l,0) + /(0,l)], 
where Wk,o = 1 - Ufc.o - Vkfi- 
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Figure 1 . Nodes of Gauss-Lobatto cubature rules of degree 5 



Both examples are constructed by following the steps in the algorithm. The 
symmetry makes the construction much easier, since we only need to consider poly- 
nomials that are symmetric under the symmetric group of A. In particular, the 
cubature rule for Wi_i^i{x, y) in Step 1 can be found by solving the reduced moment 
equations of symmetric polynomials. We shall skip details and only list the nodes 



and weights of these two cubature rules as formulated in (2.14). Their nodes are 
depicted in Figure 2. 

Example 2. Symmetric Gauss-Lobatto cubature rules of degree 5 with 12 nodes. 
This formula is in the form of (2.141 with Mg = 1 and Mi = 2. The nodes and 
weights are given below: 



"1,0 



"1,0 

"1,1 

Bi = B2 = 



Jj(7-^/7), Ai = ^(14-77), 



^(21 



2l(4%/7-7' 



720' 
-"1,2 



^(21 



2l(4%/7-7' 



720 



7 + 4V7 



'^ ~ 720' 



-V7 



D 



Example 3. Symmetric Gauss-Lobatto cubature rules of degree 7 with 18 nodes. 
This formula is in the form of (2.141 with A/q = 2 and Afi = 3. The nodes and 
weights are given below: 



"1,0 = 



"1,0 = 1^(5- V7), 
(1141-94^7); 



1 

17640 



"1.1 



= i(3-^/3 



"2,1 = 



"2,0 = "2,0 = 1^(5 + \/7), 
^2-1^(1141 + 94^7), 
\, 7.3,1 = 1(3 + 73), 



B^ = 



280' 



B2 = 



315' 



B.= 



280' 



r - 1 

•-^ - 3l5- 



D 

These cubature rules appear to be new (see the list in [1]). Their numbers of 

nodes are more than the minimal given in (2.2). The existence of higher order 



Gauss-Lobatto rules depends on the existence of minimal cubature rules of degree 
2n — 4 for the weight Wi.1^1. The latter cubature rules most likely do not exist for 
Ti > 6, but the algorithm can still be applied to produce cubature formulas with 
n — 1 points in each side of the triangle. 
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Figure 2. Nodes of symmetric Gauss-Lobatto cubature rules of 
degree 5 and 7 

It should be mentioned that there are cubature rules of degree 2n — 1 in the 



form of (2.3) that have fewer than n — 1 points on each side, naturally with more 



interior points according to (2.5), such rules cannot be constructed directly by our 



algorithm. It is possible to modify the algorithm, however, since Ci for such rules 
cannot be positive definite on n„_i but it must be positive definite on a subspace 

of n„_i. 

Finally, let us mention that our algorithm can also be modified for constructing 
cubature rules of degree 2n. This means a cubature rule of degree 2n — 3 for 
VFq+i, ^+1,7+1 in Step 1, and a quadrature for Li of degree 2n — 2 in Step 2. The 
quadrature of degree 2n — 2 is generated by a quasi-orthogonal polynomial of the 
form g„ := pn+otpn-i, where a is a free parameter which can be fixed by requiring, 
say, 5„(l/2) = 0, which means fixing the middle point on the corresponding side 
of the triangle as a node of the cubature. We have tried this construction for 
cubature rules of degree 6 with 4 interior points, 3 on each sides and 1 at each 
vertex, starting with a cubature rule for VFi,i,i of degree 3. The Lobatto type 
cubature rule of degree 6 that we obtained, however, has three negative weights. 

Acknowledgment. The author thanks both referees for their careful reading of 
the manuscript, especially David Day for his extensive and thoughtful suggestions. 
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